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We obtain an efficient description for the dynamics of nonspinning compact binaries moving in 
inspiralling eccentric orbits to implement the phasing of gravitational waves from such binaries at 
the 3.5 post-Newtonian (PN) order. Our computation heavily depends on the phasing formalism, 
presented in [T. Damour, A. Gopakumar, and B. R. Iyer, Phys. Rev. D 70, 064028 (2004)], and 
the 3PN accurate generalized quasi-Keplerian parametric solution to the conservative dynamics 
of nonspinning compact binaries moving in eccentric orbits, available in [R.-M. Memmesheimer, 
A. Gopakumar, and G. Schafer, Phys. Rev. D 70, 104011 (2004)]. The gravitational- wave (GW) 
polarizations h+ and hx with 3.5PN accurate phasing should be useful for the earth-based GW 
interferometers, current and advanced, if they plan to search for gravitational waves from inspiralling 
eccentric binaries. Our results will be required to do astrophysics with the proposed space-based 
GW interferometers like LISA, BBO, and DECIGO. 
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I. INTRODUCTION 

Inspiralling compact binaries of arbitrary mass ratio 
moving in quasi- circular orbits are the most plausible 
sources of gravitational radiation for the first generation 
ground-based interferometric detectors [1]. The avail- 
ability of highly accurate general relativistic theoretical 
waveforms required to extract the -wesk GW signals from 
the noise-dominated interferometric data is the main rea- 
son for the above understanding. The dynamics of long 
lived and isolated compact binaries can be modelled ac- 
curately in the PN approximation to general relativity as 
point particles moving in quasi-circular orbits. The PN 
approximation allows one to express the equations of mo- 
tion of a compact binary as corrections to the Newtonian 
equations of motion in powers of (w/c)^ '-^ GM/{c?R), 
where v, M, and R are the characteristic orbital veloc- 
ity, the total mass, and the typical orbital separation of 
the binary, respectively. Recently, the orbital evolution 
of nonspinning compact binaries in quasi-circular orbits, 
under the action of general relativity, was computed up 
to the 3.5PN order in Ref. [2]. The amplitude corrections 
to the GW polarizations /i+ and are also available to 
the 2.5PN order [3]. 

However, a recent surge in astrophysically motivated 
investigations indicates that compact binaries of arbi- 
trary mass ratio moving in inspiralling eccentric orbits 
are also plausible sources of gravitational radiation even 
for the ground-based GW interferometers. One of the 
earliest scenarios involves Kozai oscillations, associated 
with hierarchical triplets that may be present in globular 
clusters [4-7]. Last year, it was pointed out that dur- 
ing the late stages of black hole-neutron star (BH-NS) 
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inspiral the binary can become eccentric [8]. This is be- 
cause in general the neutron star is not disrupted at the 
first phase of mass transfer and what remains of the neu- 
tron star is left on a wider eccentric orbit from where it 
again inspirals back to the black hole. This scenario was 
very recently invoked to explain the light curve of the 
short gamma-ray burst GRB 050911 [9]. Another sce- 
nario, reported in Nature, suggests that at least partly 
short GRBs are produced by the merger of NS-NS bi- 
naries, formed in globular clusters by exchange interac- 
tions involving compact objects [10]. A distinct feature 
of such binaries is that they have high eccentricities at 
short orbital separation [see Fig. 2 in Ref. [10]]. Com- 
pact binaries that merge with some residual eccentrici- 
ties may be present in galaxies too. Chaurasia and Bailes 
demonstrated that a natural consequence of an asymmet- 
ric kick imparted to neutron stars at birth is that the ma- 
jority of NS-NS binaries should possess highly eccentric 
orbits [11]. Further, the observed deficit of highly eccen- 
tric short-period binary pulsars was attributed to selec- 
tion effects in pulsar surveys. The authors also pointed 
out that their conclusions are applicable to BH-NS and 
BH-BH binaries. Yet another scenario that can create 
inspiralling eccentric binaries with short periods involves 
compact star clusters. It was noted that the interplay be- 
tween GW-induced dissipation and stellar scattering in 
the presence of an intermediate-mass black hole can cre- 
ate short-period highly eccentric binaries [12]. Finally, a 
very recent attempt to model realistically compact clus- 
ters that are likely to be present in galactic centers indi- 
cates that compact binaries usually merge with eccentric- 
ities [13]. These above mentioned scenarios force us to 
claim that compact binaries in inspiralling eccentric or- 
bits are plausible sources of gravitatinal waves even for 
the ground-based GW interferometers. 

In order to do astrophysics with the proposed space- 
based GW interferometers, LISA [14], BBO [15], and DE- 
CIGO [16], it is required to have highly accurate GW 
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polarizations, /i+ and hx, from compact binaries of arbi- 
trary mass ratio moving in inspiralling eccentric orbits. 
Recall, the earlier discussions also indicate that stellar- 
mass compact binaries in eccentric orbits are excellent 
sources for LISA. Furthermore, it is expected that LISA 
will "hear" gravitational waves from intermediate-mass 
black holes moving in highly eccentric orbits [17-19]. Fi- 
nally, several papers which appeared recently in the arXiv 
indicate that supermassive black-hole binaries, formed 
from galactic mergers, may coalesce with orbital eccen- 
tricity [20-24]. It is interesting to note that these inves- 
tigations employ different techniques and astrophysical 
scenarios to reach the above conlusion. 

The above mentioned astrophysically inspired inves- 
tigations motivated us to extend the phasing formal- 
ism, developed and implemented with 2.5PN accuracy in 
Ref. [25], to the next FN order, namely, the 3.5PN order. 
The phasing formalism provides a method to construct, 
almost analytically, templates for compact binaries of ar- 
bitrary mass ratio moving in inspiralling eccentric orbits. 

We recall that accurate templates for the detection 
of gravitational waves require phasing, i.e., an accurate 
mathematical modelling of the continous time evolution 
of the GW polarizations. In the case of inspiralling eccen- 
tric binaries, the above modelling requires the combina- 
tion of three different time scales present in the dynam- 
ics, namely, those associated with the radial motion (or- 
bital period), advance of periastron, and radiation reac- 
tion, without treating radiation reaction in an adiabatic 
mannner. In Ref. [25], an improved method of variation 
of constants was presented to combine these three time 
scales and to obtain the phasing at the 2.5PN order. We 
note that the techniques adapted in Ref. [25] were in- 
fluenced by the mathematical formulation, developed by 
Damour [26-28], that gave the (heavily employed) accu- 
rate relativistic timing formula for binary pulsars [29, 30]. 

It is possible to extend the phasing to the 3.5PN or- 
der, mainly because of the recent determination of the 
3FN accurate generalized quasi-Keplerian parametric so- 
lution to the conservative dynamics of nonspinning com- 
pact binaries of arbitrary mass ratio moving in eccentric 
orbits [31]. This parametrization, presented in Ref. [31], 
allows one to solve analytically the 3PN accurate conser- 
vative dynamics of nonspinning compact binaries, com- 
puted both in ADM-type coordinates [32] and in har- 
monic coordinates [33], indicating the deterministic na- 
ture of the underlying dynamics [34]. Further, we recall 
that Ref. [31] extends the quasi-Keplerian parametriza- 
tion developed by Damour and his collaborators [35-37], 
which is crucial to construct the timing formula relevant 
for relativistic binary pulsars [29, 30]. This observation 
clearly reveals the Hnk, as envisaged by Damour, connect- 
ing GW observations of inspiraUing compact binaries to 
the timing of binary pulsars. In Ref. [25], explicit com- 
putations to realize the phasing at the 2.5PN order were 
done in ADM coordinates as at that time the 2PN ac- 
curate generalized quasi-Keplerian parametrization, re- 
quired by the method of variation of constants, was only 



available in ADM coordinates. However, in this paper, 
computations required for the phasing at the 3.5PN order 
are done in harmonic coordinates. Harmonic coordinates 
are preferred as computations that lead to ready to use 
search templates for compact binaries in quasi-circular 
orbits are usually performed in harmonic gauge [2]. 

This paper has the following plan. In Sec. II, we out- 
line the procedure, detailed in Ref. [25], to perform the 
phasing. Section III provides explicit formulae required 
for the 3.5PN accurate phasing, which also include 3.5PN 
accurate equations for the secular and periodic variations 
of the orbital elements involved in the phasing. The pic- 
torial representation of the main results and the related 
discussions are presented in Sec. IV. Finally, in Sec. V, 
we give a brief summary and point out possible exten- 
sions. Appendices A and B deal with computational de- 
tails and some supplementary results. 

II. BRIEF DESCRIPTION OF GW PHASING 
FOR INSPIRALLING ECCENTRIC BINARIES 

In this section, we summarize the basic ideas of GW 
phasing, as detailed in Ref. [25]. We also discuss what is 
required for this purpose and how the method of variation 
of constants will be invoked. 



A. Basic structure of GW phasing 

Let us first describe the method, developed in Ref. [25], 
to implement the GW phasing for compact binaries in 
inspiraUing eccentric orbits. The theoretical templates 
required by the GW interferometers consist of the two 
independent GW polarization states /i+ and ■ The ap- 
propriate expressions for /i+ and hx, expressed in terms 
of the binary's intrinsic dynamical variables and location, 
are given by 

^+ = ^{PiPj - mj)hj^ , (la) 

= ^{PiQ] +P]qi)hJ^ , (lb) 

where hj^"^ is the transverse-traceless (TT) part of the 
radiation field. The two orthogonal unit vectors p and 
q span the plane of the sky, i.e., the plane transverse to 
the radial direction linking the source to the observer. 

The TT radiation field is given, by the existing GW 
generation formalisms [2, 3, 38-40], as a PN expansion in 
(v/c). In this paper, for simpHcity, we will restrict hj^^ to 
its leading "quadrupolar" order and denote it by hj^^lq. 
However, higher-PN corrections to hj^^ are available in 
the existing literature [40, 41]. The explicit expression 
for hjj^lq, in terms of the relative separation vector r 
and the relative velocity vector v, reads 

,TTi 4GAi_ ,,,,, / GM \ 

hmlQ'^ -^^^^Jk,n{N) \ ll,Vj ^'^i^jl , (2) 
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where Vijkm{N) is the usual TT projection operator pro- 
jecting normal to TV, where N = R' / R' is the line-of- 
sight unit vector from the binary to the observer, and 
R' = \R'\ is the corresponding radial distance. The re- 
duced mass of the binary ^ is given by ■mim2/M, where 
M = mi +m,2 is the total mass of the binary consisting 
of individual masses mi and TO2. The components of the 
unit relative separation vector n = r/r, where r = \r\, 
and the velocity vector v — dr/dt are denoted by rii and 
Vi, respectively. 

In order to compute /i+|q and /ix|q, the expressions 
for the GW polarization states when their amplitudes are 
restricted to the leading quadrupolar order, one needs to 
choose a convention for the direction and orientation of 
the orbital plane with respect to the plane of sky. We fol- 
low the convention used in Refs. [25, 39, 40] for choosing 
the orthonormal triad (p, q, N) — N from the source to 
the observer and p toward the correspondingly defined 
ascending node — and use r = r cos (jip + r sin 0(cos iq + 
siniTV), where i denotes the inclination angle of the or- 
bital plane with respect to the plane of the sky. This 
leads to the following lowest-order contributions to the 
two independent GW polarization states as functions of 
the relative separation r and the true anomaly (j), i.e., 
the polar angle of r, and their time derivatives f and (j), 
which read [25] 



-^|(1 + C2) 2fr^sin20 

GM 9 • 9 . 9 \ „ , 

h - r cos2(/) 



'GM 



2 12 -2 

r m — r 



(3a) 



2GnC 



c-^R' 



+ r^0^ - I sin 20 



2fr0 cos 2(f> 



(3b) 



r(t), and 4>{t), relies on the possibility of going to a suit- 
able defined center-of-mass (COM) frame. Though such 
a frame exists only to the 3PN order [42], we will not 
worry about the associated "recoil" of the center of mass 
at the 3.5PN order. This is because, as demonstrated 
in Ref. [25], the influence of the recoil on waveforms ap- 
pears at the 4PN order and we may safely neglect it for 
the present computation. 

A second point, again detailed in Ref. [25], is that the 
GW polarizations /i+lq and /ix|q only depend on the 
temporal evolution ofr{t), (t>{t)^ f(i), and (/)(i), as we are 
dealing with nonspinning compact objects. In the pres- 
ence of spin interactions, the orbital plane is no longer 
fixed in space and one needs to introduce further vari- 
ables [30, 36, 43, 44]. However, expressions for /i+|q and 
/ixIq, valid for spinning compact binaries are recently 
obtained in Ref. [43] and the phasing for such binaries 
will be reported elsewhere. 

Finally, we want to discuss the choice of the under- 
lying coordinate system. Though the expHcit functional 
forms of /),+ (r, 0, f, 0)|q and ft-x (r, 0, f, 0)|q, as well as 
the phasing relations r(t), (j){t), r{t), and (j){t) depend on 
the used coordinate system, the final results h+{t) and 
hx{t) do not. Note that hj^"^ and therefore h^{t) and 
hxit) are coordinate independent asymptotic quantities. 
In this paper, we consistently work in harmonic coordi- 
nates because (i) they will allow us to write down explicit 
analytical expressions for the orbital phasing r{t), 4'{t), 
r{t), and (/)(i), based on Ref. [31], and (ii) the harmonic 
coordinate systems are used in the standard GW gen- 
eration formalisms to derive the amplitude expressions, 
giving higher-PN corrections to hj/j^ . 

In the next subsection, we will explain a version of the 
general Lagrange method of variation of arbitrary con- 
stants, which was employed to compute, within general 
relativity, the orbital evolution of the Hulse- Taylor bi- 
nary pulsar [27, 28], and appHed for the 2.5PN accurate 
phasing in Ref. [25]. 



where C and S are shorthand notations for cos i and sin i, 
respectively. The orbital phase is denoted by 0, = 
d(j)/dt, and r = dr/dt = n ■ v. 

In order to achieve the GW phasing for compact bi- 
naries of arbitrary mass ratio moving in inspiralling ec- 
centric binaries, we need to provide explicit expressions 
describing the temporal evolution of r{t), (j){t), r{t), and 
4>{t), the only dynamical variables appearing in Eqs. (3). 
Following Ref. [25], we refer to, as phasing, an explicit 
way to define the latter time dependences. This is the 
crucial input to derive ready to use waveforms h^{t) and 
hx(t). Adapting the formalism presented in Ref. [25], 
we will provide the above desired time evolution to the 
3.5PN order in an almost parametric manner. 

However, there are a few points that we want to clar- 
ify before we describe the above procedure. First, as ex- 
plained in Ref. [25], the possibility of obtaining explicit 
expressions for the GW polarizations /i+lq and h^lq in 
terms of the relative dynamics, specified by r{t), 4'{t), 



B. Improved method of variation of constants and 
its implementation at the 3.5PN order 

The method, presented in Ref. [25], begins by split- 
ting the relative acceleration A. of the compact binary 
into two parts, an integrable leading part Aq and a per- 
turbative part A! , as A = Aq + A! . The method first 
constructs the solution to the "unperturbed" system, de- 
fined by 



r = V , 

i) = Ai{r,v). 
The solution to the exact system 



i} — A{r, v) . 



(4a) 
(4b) 



(5a) 
(5b) 
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is then obtained by varying the constants in the generic 
solution of the unperturbed system, given by Eqs. (4). 

In this paper, we work to the 3.5PN order and therefore 
Aq will be the acceleration at the 3PN order and A! con- 
tains the reactive 2.5PN and 3.5PN contributions. The 
method assumes — as is true for — that 

the unperturbed system admits sufficiently many inte- 
grals of motion to be integrable. At the 3PN order with 
Ao = ^3PN, we have four first integrals: the 3PN ac- 
curate energy and the 3PN accurate angular-momentum 
vector of the binary. We denote these quantities, written 
in the 3PN accurate COM frame by ci and Cj: 



ci =£^(r,t;)|3pN COM' (6a) 
4 =^'(^''")l3PN COM- (6b) 

At the 3PN order, the functional form of the solution to 
the unperturbed equations of motion in the COM frame, 
following Refs. [27, 31], may be expressed as 



r = S{1; ci, C2) , 
dS 

r = n—[l;ci,C2) , 

(/)= X + W{l;ci,C2), 
dW 

(j)^ {1 + k)n + n—-{l; Ci, C2) , 
ol 



(7a) 
(7b) 
(7c) 
(7d) 



where A and / are two basic angles, which are 27r periodic, 
and C2 = |c||. The functions S{1) and W{1) and hence 
their derivatives ^(0 and ^j^(0 are also periodic in 
I. In the above equations, n denotes the unperturbed 
"mean motion", given by n = 27r/P, P being the radial 
(periastron to periastron) period, while k = A$/(27r), 
where A<i> represents the advance of periastron in the 
time interval P. The explicit 3PN accurate expressions 
for P and k in terms of ci and C2 are obtainable from 
Refs. [31, 45]. The angles I and A satisfy, still for the 
unperturbed system, I = n and A = (1 -|- k)n, which 
integrate to 



I = n{t- to) + Q , 
A = (1 + k)n{t - to) 



(8aj 
(8b) 



where is some initial instant and the constants c; and 
c\ are the corresponding values for I and A. Note that the 
unperturbed solution depends on four integration con- 
stants: ci, C2, c/, and c\. 

Before we explain the prescription to obtain the phas- 
ing at the 3.5PN order, let us first present our description 
of the 3PN accurate conservative dynamics. The 3PN ac- 
curate parametric solution to the dynamics of compact 
binaries of arbitrary mass ratio moving in eccentric or- 
bits, derived in Ref. [31], allows us to describe the orbital 
motion in the 3PN accurate COM frame. The vectorial 
structure of [see Eq. (6b)], indicates that the unper- 
turbed 3PN accurate motion takes place in a plane, the 
so-called orbital plane. The problem is restricted to a 
plane even in the presence of radiation reaction. There- 
fore, we introduce polar coordinates in the orbital plane. 



r and 4>, such that the relative separation vector r takes 
the following form r = r cos <j)i + r sin (f>j , where we may 
choose i = p and j = cos iq + sin iN as the definition 
for the basic vectors. The dynamical variables r and ip, 
available in Ref. [31], are parametrically given by 



(1 — e, 
(1 + k)v 



(9a) 



./40 



^4 



h<t> 
sin 3w - 



sin 2v 



sm 4i; 



sin 5u , 



where v = 2 arctan 



1 



1/2 



tan ■ 



(9b) 
(9c) 



In the above equations, a^, e^, and are some 3PN 
accurate semimajor axis, radial eccentricity, and angu- 
lar eccentricity, respectively, while /40, /e^, g4</-, 560, «6<^, 
and are some PN accurate orbital functions. The 
following 3PN accurate Kepler equation connects the ec- 
centric anomaly u (and hence the true anomaly v) to the 
coordinate time t and reads 



I = n {t — to) 



u — et sin u + 
fit , fet 



+ ) (^' - ") 



s- sm 6v . 

c6 



«6t 



sin 2v 



(10) 



Here et denotes some time eccentricity and 1744, get, fit, 
fet, iet; and h^t are some PN accurate orbital functions. 
All these PN accurate orbital elements and functions, 
expressible in terms of ci, C2, and 77 = /i/M, are available 
in Ref. [31]. 

With the help of the above 3PN accurate parametric 
solution, we write down expHcit expressions for the func- 
tions S{1) and W{1): 

S {I; d, C2) ^ ar {1 — er cos u) , (Ha) 
W{1; ci,C2) ^ {1 + k)iv - I) + I H + ^ \ sin 2v 



.940 
,4 



c 



.960 



sinSw 



■ sin 4w 



sin 5u . 



(lib) 



We emphasize that for the phasing the anomalies v and 
u in Eqs. (11) have to be expressed as functions of Z, ci, 
and C2 . This can be achieved in the following way, written 
symbolically as u = V(/;ci,C2) — V[U{l;ci,C2)] and u = 
U{1; ci, C2). First, recall that the function v = V{u) is 
defined by Eq. (9c), and second, the function u =U{1) is 
defined by inverting the Kepler equation / = /(w), given 
by Eq. (10). Finally, the function v = V{1) is obtained by 
inserting u = U{1) in w = V{u), i.e., v = V{1) = V[U{1)]. 
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In our computations, we use the following exact rela- 
tion for V — u, which is also periodic in u: 



2 tan 



1 



, cos u 



(12) 



where 



(1 



1 



e|)/e0 [for the derivation of 

Eq. (12), see Appendix A]. 

In line with the method of variation of constants, we 
write down the following general solution to the reactive 
3.5PN accurate dynamics, Eqs. (5), as 

r = S{l;ci,C2) , (13a) 

r = n—{l;ci,C2) , (13b) 

= A + Ty(/;ci,C2), (13c) 
dW 

(/)= {l + k)n + n—{l;ci,C2). (13d) 

However, the constants ci and C2, appearing in Eqs. (7), 
are now functions of time in Eqs. (13): ci ~ ci(0 and 
C2 = C2 {t) . The temporal variation of the basic angles I 
and A, entering Eqs. (13), is now given by 



I 



ndt + ci{t) , 



A= / {I + k)ndt + cx{t) , 

J to 



(14a) 
(14b) 



involving two new temporally evolving quantities ci{t) 
and cx{t). In the method of variation of constants, we 
search for solutions of the exact system, Eqs. (5), in the 
form given by Eqs. (13) and (14) with four "varying con- 
stants": ci(t), C2(i), Q(i), and CA(i). These four new vari- 
ables replace the original four dynamical variables r{t), 
r{t), 4>{t), and 0(t), and satisfy, Hke the original phase- 
space variables, first-order evolution equations [27, 28]. 
These first-order evolution equations for Ci (t) , C2 (t ) , ci{t), 
and c\{t) are available in Refs. [25, 28] and read in our 
notation 



dci 
~dt 
dc2 
~dt 
dci 
dt 
dcx 
dt 



dci{r,v) 
dc2{r,v) 



'dl 

dWdci 



A'- 



dS dci 
dci dt 
dW dci 



dS dc2 
dc2 dt 
dW dc2 



dl dt dci dt dc2 dt 



(15a) 
(15b) 
(15c) 
(15d) 



In addition, an alternative expression for dci/dt is pre- 
sented in Refs. [25, 28], which we used for some supple- 
mentary checks. 

We observe that the definition of the sole angle I, ap- 
pearing in Eqs. (15), given by I — J^^ n[ca{t)]dt + ci{t), 
where c^, a = 1, 2 stands for ci and C2, is equiva- 
lent to the differential form dl/dt — n{ca) + dci/dt — 



n'(ca) + Fi{l;Ca). This allowed Ref. [25] to write the above 
set of equations, namely, Eqs. (15) for Cq as functions of 
t, as a set of differential equations for Cq as functions of 
which symbolically reads 



dCa FQ^{l^Ca) 

dl n{ca) + Fi{l; Ca) ' 



(16) 



where a = 1, 2, Z, A and a = 1, 2. Neglecting terms 
quadratic in F^, i.e., quadratic in the perturbation A!, 
terms of 0{c~^^) (and higher-PN orders), we can simplify 
the system above to 



dca 

~dr 



■n{ca) ~ 



Ga(l] Co) 



(17) 



It is important to note at this stage that in the above pre- 
scription, which neglects 0(0"^°) terms, the right-hand 
side of Eq. (17), first, is just a function of ci, C2, and the 
sole angle I (and not of A), and second, it is a periodic 
function of I. This periodicity, together with the slow 
evolution of the Cq's, clearly indicates that the evolution 
of Cq(/) contains not only a slow secular drift, but also 
fast periodic oscillations. The secular drift occurs in the 
radiation-reaction time scale, while the time scale for the 
periodic oscillations is that of the orbital motion. 

Following Ref. [25] , we model the combination of the 
slow drift and the fast oscillations in Cq(Z), by a two-scale 
decomposition of Ca{l), which reads 



,(0 = Ca(0 + Cq(0 ■ 



(18) 



In the above equation, Ca{l) denotes the slow drift, which 
accumulates over the radiation-reaction time scale to in- 
duce large changes in Cq(/), while the fast oscillations in 
Ca{l) are denoted by Ca{l), which will be always smaller 
than Ca{l). 

In this paper, we are interested in the 3.5PN accurate 
evolution of Ca{l) and Ca{l)- The corresponding evolution 
equations for Ca(l) and Ca{l) follow from the combination 
of Eqs. (41), (43), and (44) of Ref. [25] and read 



dCa 

~dr 

dCa 



G^ica)^^ f ^ G{l;ca)dl, (19a) 

Jo 



= G^il;ca) ^ G^{l;ca) - G^iCa) . (19b) 

Finally, the arguments that gave the 2.5PN accurate 
unique zero-average expressions for Cq-(Z) in Ref. [25], can 
be extended to obtain the 3.5PN accurate expressions for 
Ca(/). This leads to the solution of Eq. (19b), considered 
for fixed values of Ca , 



Ca{l) 



Ga{l;ca)dl 



F^{l;ca)-. (20) 
n 



We recall from Ref. [25] that the indefinite integral in 
Eq. (20) is defined as the unique zero-average periodic 
primitive of the zero-average (periodic) function Gail)- 
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During that integration, the arguments Ca are kept fixed, 
and, after the integration, they are replaced by the slowly 
drifting solution of Eq. (19a). 

Furthermore, a consistency check of the two-scale 
method was also provided in Ref. [25] by exploring the 
effects of the above neglected second-order terms in 
Eqs. (17)-(20). It was shown there that the separation 
between the two scales remains valid on very long time 
scales and that second-order (and higher-order) effects 
cause only fractionally small separate corrections to the 
evolution of Ca{l) and Cq(/). 

In the next section, we will implement the arguments 
developed above and obtain explicit expressions to per- 
form the GW phasing at the 3.5PN order. 

III. ANALYTICAL EXPRESSIONS FOR THE 
3.5PN ACCURATE PHASING 

In this section, we apply the above described improved 
method of variation of arbitrary constants, which gave 
us the evolution equations for Cq, and to achieve the 
GW phasing in the following way. First, we compute the 
3PN accurate parametric expressions for the dynamical 
variables r, r, (j), and (/> entering the expressions for /i+|q 
and hxlq, given by Eqs. (3). Then we solve the evolution 
equations for ci, C2, Q, and c\, given by Eqs. (17), on the 
3PN accurate orbital dynamics, given in Eqs. (7). This 
leads to an evolution system, given by Eqs. (19), where 
the right-hand side contains dominant 0{c~^) terms and 
their first corrections, i.e., 0{c~^) terms. Later, we will 
impose these variations on the 3PN accurate expressions 
for the dynamical variables r, r, (j), and 4>, appearing in 
/i+Iq and hxlq, given by Eqs. (3). This allows us to 
obtain GW polarizations, which are Newtonian accurate 
in their amplitudes and 3.5PN accurate in the orbital 
dynamics. Following Ref. [25], the above procedure is 
called 3.5PN accurate phasing of gravitational waves. In 
our computations, the quasi-periodic oscillations in Cq, 
governed by Ga, are restricted to the IPN reactive or- 
der. Therefore, we will not explore higher-PN correc- 



tions to the above Ga- However, in Appendix B, we 
will present the consequences of considering higher-PN 
corrections to Ga by computing 0(c~^) contributions to 
relevant dca/dt. This is desirable as Ga directly con- 
tributes to the highly important adiabatic evolution of 
/i+ and hx ■ 

In this paper, we follow Ref. [25] and employ as ci the 
mean motion n, and as C2 the time eccentricity et, instead 
of the energy E and the angular momentum L, respec- 
tively, to describe the PN accurate reactive dynamics. 
This can be done by employing 3PN accurate expressions 
for n and et in terms of E and L, derived in Ref. [31]. This 
implies that first, we have to express the 3PN accurate 
orbital dynamics in terms of I, n, and et- Second, using n 
and et , instead of E and L as ci and C2 , respectively, we 
need to derive the evolution equations for dn/dt, det/dt, 
dci/dt, and dc\/dt in terms of ^, n, and e*. This will 
follow straightforwardly from Eqs. (15). Using these ex- 
pressions, the evolution equations, namely, Eqs. (19), for 
n, et, ci, Ca, n, et, h, and c\ will be obtained in terms of 
I, n, and et- 

Let us now obtain the explicit 3PN accurate expres- 
sions for the orbital dynamics. 



A. 3PN accurate conservative dynamics 

As mentioned earlier, we restrict in this paper the con- 
servative dynamics to the 3PN order. Below, we present 
the 3PN accurate orbital dynamics, namely, r, f, 4>, and 
(j), as given by Eqs. (13), explicitly in terms of u, n, and 
et- This straightforward computation employs explicit 
expressions for the PN orbital elements a^, er, e^, k, f^^, 
fsrji, 540, 560, isip, hQ^, g^t, get, fit, fet, iet, and het of 
the generalized quasi-Keplerian representation, available 
in Ref. [31]. We note that the above PN orbital elements, 
given in terms of E and L in Ref. [31], can easily be ex- 
pressed in terms of n and et with the help of the following 
3PN accurate relations for —2E and —2EL'^: 



-2EL^ 



2E = (GMnf^H 1 + ^(15 - 77) + V 



^2/3 



^4/3 



12 



24 



- 3577^ 



864 



= (15 + 2377-2077^ 



15 — 15?7 — 77 
18 

^ (1 - ^2)3/2 



24 



= (5-2r,) 



5184 



4995 - 6075?7 - 45077^ 



(11520 - 1596877 + 1237r^77 + 201677^ 



(21a) 



(1 



2/3 



4(1 -e?) 



[9 + 77 - (17 - 777)6?] 



4/3 



24(1 -e?)2 



2\„2 
t 



(225 - 277?7 + 2977^)6^ - (360 - 14477)e^/l - e: 



189 - 45?7 + 77^ - 2(111 + 777 + 1577^)e 

e 



6720(1 -e?)3 



35(5535 - 9061?7-f 2467r^77 



14277^ - 77-^) + (299145 - 119766777 + 258307r^77 + 17325077^ + 2M5r]^)ef + 35(3549 - 1278377 -I- 615477 
13177^)6^ - 35(2271 - 73817/ + 241477^ - 65jf)e^ + 70 [24(45 - I377 - 27;^) - (17880 - 2OI2O77 



7 



^UStt^j] + 22h&ff)e1 + 96(55 - 40/; + ■iff)ef\ y'l -e^j^ , 



(21b) 



where ^ = GMnjt? and ?/ = [ijM = mYmij (m\ + TO2)^- These two relations follow from inverting the 3PN accurate 
relations for the orbital period P = 27r/n and the squared time eccentricity in terms of E and L presented in 
Eqs. (25c) and (25d) in Ref. [31]. 

In addition, to compute expressions for r and 0, we use the following relations: 



a„e„ smu 



dl ' 



dS_ . du 

'dl 
dW 

'df 

du 
'dl 

dv (1-6^)1/2 



1 + ^. + 21":^ + /^Vos2z; + 3(^^ + cos3« + 4^ cos4z; + 5^ cos5^ 



1 — ct cos u 



54t get 



(22a) 

dv du ,^ , . , , , 
---(1 + .), (22b) 



54t 










C6 







COS V + 2— COS 2v + 3—2- cos 3v 



dv 
du 



(22c) 
(22d) 



1 

The radial motion, defined by r{l, n, et) and f{l, n, et), reads (both in the compact and in the 3PN expanded form) 

r = S{l,n,et) = ar(n, et)[l - Crin, et) cosu] = tn + riPN + r2PN + t'spn , (23a) 

where 

1/3 



[ ) ^ etcosii) , 

£2/3 

fiPN = 7'N X 777:; ^[^18 + 277 — (6 — 777)64 cosw] , 



f 2PN = X 



(1 — et cos u) 

£4/3 



72(1 - e2)(l - et cosw) 
36(5 - 2?/) (2 + et cosu)^/! - e^j , 

^2 / 70 56221 123 



I -72(4 - 77]) + [72 + 30?7 + Sr;^ - (72 - 231?; + 35?72)et cosw] (1 - e^) 
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(23c) 

(23d) 
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(23e) 



and 



9S n ^ • • 

r = n—{l,n,et) = 7-n + tipn + r2PN + rsPN ■ 
dl 



(24a) 



where 
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(24e) 



The angular motion, described in terms of (j> and (p, is given by 

<^(A,0 = X + W{1), 
A = (1 + k)l , 
W{1) = + WiPN + VF2PN 



3PN 



(25a) 
(25b) 

(25c) 



where 



W2PN 
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1 



4(1 



o2\2 



[78 - 2877 + (51 - 26?7)e?] 



128(1 -e?)3 



1 18240 - 2537677 + 4927r277 + 8967?^ 



(28128 - 2784077 + 12377^77 + 5120772)e? + (2496 - I76O77 + 1040772)e^ 



+ [1920 - 76877 + (3840 - 153677)e?] y^l - e^j 

Tj — + et sin u , 
3^2/3 



1 



e^ 
^t 



(v — u + et sin u) . 
^4/3 



78 - 2877 + (51 - 2677)e2 - 6(5 - 277)(1 - e^)^/^] (v ~ u){l - et cos 



32(1 - e?)2(l - et cos7i)3 

+ |624 - 28477 + 477^ + (408 - 8877 - 8772)e^ - (6O77 - 477^)6^ + [-1872 + 79277 - 877^ - (1224 - 38477 
-1677^)64 + (I2O77- 8772)et] etCOSu+ [l872 - 73277 + 477^ + (1224 - 50477 - 8772)e2 - (6077-477^)64] 
X (et cosu)^ + [-624 + 22477 - (408 - 20877)6^^] (et cosu)3|et sinu + | - (8 + I5377 - 2777^)6^^ 



(25d) 
(25e) 

(25f) 



9 



(1-6^)3(1 



285 793 123 2 „ 2 
77+ —TT j] + 7ri 




- 677^)6^] et cos M + [-8 - 148?7 + Urj'^ - (?7 - Sry^). 



(25g) 



^879 435 123 , 

n H TT 77 + 40?7' 

V 4 2 ' 128 ' ' 



.2 \ „2 



''959 118663 41 , 671 , 
-^4^'^-128"^+l^'^ 



"185 



157„2 , 37^3' 



■^"64 V ^ u u. 

r265 126159 287 

u)(l — et cosM + < 77 H 

I 2 560 ' 64 

19 23967 41 . 525 

W+64"" + T6 



41 2 2 A65 451 41 2 

287 o 133 ., 



2 ^''^ 



1 3 

24 ' 



48 ' 



16 



4715 713597 41 2 2987 2 35 ,\ 2 

13 3 



525 2 11 3\ 4 

1345 75417 697 2 2 2 o 

+ -^^V - -^TT ?7 - 6577 - -77 



^ 541 




840 

173 2 23 3^ 
~^ ~ 24" J 



70 

41 2 419 2 3 3\ 4 
■TT ?/ ^?7 - l^t 



4 



4981 2 71 3 /4635 1609903 
48 " +48^ V~2 840 
1077 2 29 o\ 4 /1585 



^23 73 2 
.24^'" 24'' 



115 80433 _ 
~ + 280 16" ■' 4 

3493027 1353 , 

77 H TT 77 



et cos u 



1365- 



^ 24 



64 

145 2 2 3\ 6 
1^" " 3" I 



135 



1385 - 



4 



1701499 2 2065 2 9 3 

71 — 41tT 7) 77 77 

840 ' ' 24 ' 8 ' 

' 62723 41 2 755 27 3' 

r'^'+W" 16" 

'1405 837281 1271 . 1807 2 13 3 

77 H TT 77 H 77 H ?7 

2 840 ' 64 ' 48 ' 48 ' 

"175 20887 41 2 637 2 13 3\ 4 
-''+64"'' + T^'' +16''T- 



123 o 3881 o 25 q \ n 

■ 23 73 2 13 3 
"48^+48'' ^48^ 

287 2 



1680 
12769 



32 



35 



123 2 



280 

123 2 -7 2 
— -TT rj — irj - 



32 



15 - 677 - 



"91 



412487 615 2 
^36r''+ 256""" 



- 4077^ 



879 435 123 , 

1 77 7r"77 

42' 128 ' 

2845 2 113 3 
77 if 



4555 578789 
"^ + ^8(r''" 64 " '' 8 

521 35 , 9 3\ ^, 
^^^^ M^" ^ S''^) (etcosii)- 

4475 150683 451 . ''^ ^'^ 

I 77 H TT 77 

V 4 140 ' 128 ' 

1 83 2 13 3 \ K / 

— 77^ 77'' I et [et cosm) 



(etcosii)^ 
2955 o 27 3^ 



48 48' 



3143 2 13 3 

V V 

16 ' 16 '/ 

793 



285 



39 55 65 2 \ 4 . 
-y + -^--'7 ]et {et 



64 ^ 64 



3 173297 41 2 

1 7? TT 77 

, 2 4480 ' 64 ' 



cosw)'^|'et sinu 



7319 2 
"384"'' 



209 3^ 

395 2 33 

?7 + -77 

12 ' 4 ' 

35 3 

-16'' 



4 



115 325 2 21 3\ 6 /I 
""2?'' + Is"" ~ le" j + V 16" ^ le" 

/ 523 22553 779 2 3041 2 63 3\ . 

1 77 TT 77 77 H 77 e: 

V 4 70 ' 128 ' 32 ' 32 ' ' ' 

533 



4480 

5 2 5 3\ 8 r 275 49477 287 2 

—77'^ H 77-^ e? + ! 77 7r^77 

i« ' 16 ' y 4 420 ' 128 ' 

67 2 



5 52821 41 2 

77 H TT 77 H 77^ 

2 1120 ' 32 ' 12 ' 



^53 



1127 204563 41 2 161J 
56^"+¥""+^ 



et cos u 



121 2 3 . 

^" - 8" 

41 2 1615 21 3' 



75337 779 2 3269 2 5 3 
W+T28"" + ^" +16" 
14869 41 2 1443 2 147 3' 
^"-64"" + -6r" -"6?" 



23 



1 9 



16 ' 



1169 52937 205 . 

1 77 TT 77 

4 210 ' 128 ' 



(et cosm)^ - 

2291 2 



541 18449 697 . 

1 77 7r^77 

4 84 ' 128 ' 24 



sT" +^"V'*'"l2 + ^ 

221 2 33 /297 110277 41 , 

?7 + -77 + 77 H TT 77 

24 ' 8 ' V 2 1120 '256 ' 



4 

165 3 



2240 

697 1073 2 9 . 



^4 



\4 



128 
'1 497 



283 114761 205 , 221 2 3 3 / 

77 H TT 77 H ?7 + -77'^ + 

4 1680 '128 ' 24 ' 8 ' V 2 1120 



533 o 1 3^ 
■"-1^" +2", 



4 



{et cosu)3 



1439 2 87 

77^ 



64 " 64 



10 



128 



61 

384' 



19 
128 



{et cosm)^ — [l5 — 677 + (30 — \2ri)ef\ (e* C0SM)^|e4 sm.u^ 1 — e^^ 



(25h) 



and 



^IPN 



(26a) 



where 



^IPN 



»2PN 



(1 — et cosuY ' 



0N X 



(1 — )(1 — et cosw) 



[3-(4-ry)e? + (l-r7)et. 



(26b) 

(26c) 



^4/3 



144 - 48?/ - (162 + 68?7 - 2r^^)ei + (60 + 267/ - 2077'')e* + (18?/ + 12r;2)e'^ 



i^SPN = <^N X 



12(l-e2)2(l-et cosu)3 
+ [-216 + 12577 + Tf + (102 + 188?7 + 1677^)6? - (12 + 97?/ - Tf)ef\ et cosw + [IO8 - 97?/ - brf 
+ (66 - 13677 + ATf )ei - (48 - 17?/ + 1777^)6^] (e* cos uf + [-36 + 2?/ - 877^ - (6 - 7O77 - 1477^)6?] (e* cos uf 

+ 18(1 - et cosw)2(l - 2e1 + et cosw)(5 - 2?/)y^l -e^j , (26d) 

1447 205 2 „ 2 n 57021 205 , 361 , 1 



(1 - elY{l — et cosm)^ 



12 



205 2 „ 2 

tt"^?? + 277^ 

64 ' ' 



) el 



-50 



863 



175193 41 



840 



-77 - 
47 



-TT 77 - 



317 2 47 3 , 4 
6 ' 12 ' ' * 



18 



2987 



41 



127 



77 H TT 77 H 77" 

210 ' 32 ' 6 ' 



25 



?7 — 377 



-,10 



488539 
-285 + —TTT—V 



1025 



121 94097 



-54 
451 



2 140 
16531 



451 
"64 



4571 



■TT 77 - 



24 



73 
24' 



182 - 



840 
114683 



64 



-TT rj - 



367 



24 



205 



1987 



168 "''+1^"'^+ 24 



59 
24' 



41 



?? TT ri 

210 ' 16 ' 



769 
'24 



o 1 2 1 i-j 3 q 

TT 77 H 77 + -77-^ 

16 ' 24 ' 8 ' 



17 . 



268137 



379 25 2 



6" 



et cosM 



213 ^77 H ^7r^77 + 1^772 - -if ] 

280 ' 16 ' 24 ' 



329 2 , 45 3 



24 



54 
697 



43177 



41 



123 
— 7 
L6 

227 



7 



-?? H TT ?) ?7 77" 

840 ' 32 ' 8 ' 24 ' 



-V I ^t 
5 3 



529223 
-243+ 77 



165061 

411 77 

168 ' 

369 



3"^ 



288269 697 . 189 . 55 . , 



79717 



287 



7705 



137- 



102349 



41 



53 



55 



"42^"+¥""+24" "24^ 



18 347 



78- 



77011 
420 



-77 



451 2 281 

■TT 77 ■ 



13 



-23 



64 ' 24 ' 
5699 41 2 641 

7? TT 71 -I 77 

105 '32 ' 24 ' 



325 92555 



29 



41 2 2 2 



24 

69 13537 
y + 140 



77-^ 



2 168 
11 23 



77- 



245 2 1 ^ 
IT'' "24^^^ 



3083 



840 

119 

(et cos uj 



32 



7r^77 



(et cosm)^ 



451 2 

TT 77 H 

64 ' 24 



33 



8 



24 



29 
24' 



(et cos 7i)^ 



139 
-6 + -— ?7 



?7 



64 



■7r^77 



533 91 2 



et cosu)^ + 15840 - I6O6477 + 12377^77 + 96077^ - (38400 - 38464?; + 24677^77 + 2976?72)et 
19z I 



+ (9600 - 44I677 - 57677^)64 + [9600 - 768077 + 153677^ + (8640 - 2147277 + 24677^77 + 1344772)et 

+ (7680 - 6816?7 + 2304?;2)e^] et cos u + [-8160 + 12512?/ - 1237r2?7 - 76877^ - (4800 - 547277 + 1824772)et ] 



X (et cosu)^|(l - et cosu)^yl — e^ 



(26e) 



11 



In the above equations, the eccentric anomaly u = U{l, n, ej) is given by inverting the 3PN accurate Kepler equation, 
Eq. (10), connecting I and u. The 3PN accurate Kepler equation in terms of n and et reads 
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Evidently, Eq. (27) does not have a IPN contribution, 
similar to its ADM counterpart, given by Eq. (52) in 
Ref. [25] . Note also the expected differences in the Kepler 
equation at higher-PN orders between the harmonic and 
the ADM gauge. 

In addition to the above explicit expressions for r, r-, 
0, and (/), we also need to evaluate the right-hand side of 
Eqs. (15). The next subsection contains a sketch of these 
computations and the explicit final expressions. 

B. 3.5PN accurate reactive dynamics 

Recall that the method of Ref. [25] is general and can 
be applied, in principle, to any PN accuracy. However, in 
this paper, we study the effects of the 2.5PN and 3.5PN 
contributions to radiation reaction on the 3PN accurate 
conservative motion. Accordingly, we will truncate away 
all effects that would correspond to higher-PN orders. 



Further, it should be noted that in the case of nonspin- 
ning point masses, it will be highly difficult to go beyond 
3.5PN accuracy for the oscillatory effects associated with 
the two-scale decomposition, as given by Eq. (18). How- 
ever, secular effects can be computed to higher-PN or- 
ders, e.g., see Appendix B. 

Below, we provide the required inputs for the compu- 
tation of the 3.5PN accurate evolution equations for the 
sets {cq} and {cq.}, where a. = n, et, ci, c\. Naturally, 
we require A! to the 3.5PN order for this purpose. The 
appropriate 3.5PN accurate expression for A! has to be 
in harmonic gauge, as our conservative 3PN dynamics is 
given in the that gauge. The expression for the relative 
reactive acceleration, to the 3.5PN order, in harmonic 
gauge, available in Ref. [46], reads 

•A! = ^'2.5PN + ^'s.SPN , (28) 

where 
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where = v v = -r^ + (j)'^ . Following Ref. [25], we ob- 
serve that explicit computations of the right-hand side of 
Eqs. (15) and (17), respectively, require only IPN accu- 
rate expressions for the orbital elements. This is mainly 
because of the fact that we are trying to obtain the phas- 



ing to the 3.5PN order and the reactive dynamics only 
involves 2.5PN and 3.5PN contributions. However, this 
does not mean that the orbital dynamics is only IPN 
accurate. In all expressions where they are needed, we 
have to include the appropriate PN accurate contribu- 
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tions. The phasing formahsm allows us to impose the 
fully IPN accurate reactive dynamics on the 3PN accu- 
rate conservative dynamics to provide the 3.5PN accurate 
phasing [see Sec. V in Ref. [25] for details]. 



Finally, the evolution equations for dn/dl, det/dl, 
dci/dl, and dc\/dl in terms of u{l, n, et), n, and ej, follow 
as 



dn ( 6 32 49 - Oe^ 35(1 - e?) ) ^^/^nr? f -360 + 1176r? 2680- 1170477 



dl 5 ] 35 I X^ 

+ [ - 4012 + 34356?? + (36 - 75677)e?] ^ + [1470 - 47880?/ - (350 - 17080?7)e^] ^ + [135 1 + 3 1 7807? 

- (24220 + 305207,)e? + (10710 - 126o',)e^] \ - + 5880,)(1 - e.^)^ ^ yXim-e\f\ 

X x^ J 

det 8£5/3 (2 _ r 3 17 49 - 9e? 35(1 - e^) 1 f^^^r? f r , ^ 9n 1 

-1 = ^ — tZ J + ^ ^ „ ^ + ^ <^ - 4320 + 6636?7 + 5328 - 764477 6? — 

dl 15et X' X^ J 315et\L ' X^ 

+ [20430 - 3617677 - (25470 - 41216?7)et] ^ + [ - 73650 + IO6I2O77 + (112692 - 11620077)6^ 

- (39042 - 10080?7)6t] \ + [102312 - 15428077 - (201264 - 216160?7)6,2 + (98952 - 6188077)6^1 4r 

+ [2730 + 9534O77 + (210 - 18690077)6^ - (8610 - 87780?7)6^ + (5670 + 378077)e^] ^ 

X 

(82782 + 1764077)(l-e?)3 35280(1 -e?)^ \ 

? ? J' ^ ^ 

dci 8C^/3,/sinu [126^ 3 - 436^ -14 + 23e?-9e^ 35(1-62)2] ^7/3^ gjn ^ f (_4176 + 940877)6? 



dl 156t [ -^4 ^5 j 3156t [ X' 

+ [ - 4320 + 663677 + (23808 - 71596?7)6?] ^ + [I6IIO - 29540?? - (23862 - 128240?7)6? 

+ (4392 - 630077)6^1 ^ + [ - 57540 + 76580?? + (103320 - 159880?;)e? - (45780 - 83300?7)6fl 4r 

X X 

+ [44772 - 777OO77 - (59934 - 15162077)6? - (14448 + 70140?7)ef + (29610 - 3780?7)e^] ^ 

X 



(47502 + 17640?7)(1 - 6?)-^ 35280(1 - e?) 



X* ^9 



(30c) 



dc\ %i^l'^r\ sin u 
dl 15ef 

2^''/^?7 sin t? 



3 14 -9e? 35(1-6?) 



X"* X^ X^ 



1-ef- 



126? 3-43e? 14 - 236? + 96^ 35(1-6?) 



^3 ^4 ^5 ^6 



[1404 - 33I877 + (360 + 306677)6?] — + [ - 4527 + 14770?? - (4029 + 15232?7)6? 



3156t(l-6jj ■ "^X 

(576 + 882?;)e^] \ + [19950 - 38290?? - (38640 - 58520?;)e? + (186 90 - 20 2 307?)6^] ^ 

X X 

[ - 22386 + 38850?? + (27447 - 758107?)6? + (12264 + 35070?;)ef - (17325 - 1890?;)e^l 4 

X 



(23 751 + 88 207?)(l-e?)3 17640(1 - 6: 



?)^ 



^8 ^9 



1 - 6? + [(2448 + 47047?)6? + (2088 - 470477)6^] ^ 

X 



' - e? 



[ - 1404 + 3318?? - (4332 + 39116??)e? - (11904 - 35798??)6^] — + [4527 - 14770?? - (14190 - 78890??) 

X 

(11859 - 67 2 707?)6f - (2196 - 3150?;)6^1 ^ + [ - 19950 + 38290?? + (62790 - 118230??)6? 

X^ 

(65730 - 121590??)6^ + (22890 - 416507?)6^] ^ + [22386 - 38850?? - (52353 - 114660?;)e? 

X 
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+ (22743 - 11088077)e4 + (22029 + 33180?/)e? - (14805 - 189077)6?] ^ + (23751 + 88 2077)(1 e^) 



17640(1 - e^)^ \ 48^'^/^ri{v - u) [ 1 5_ 
? /+ 5(1 -e?) Ir^"^ 



2\4 

(30d) 



where x = 1 ~ ej cosw and u = u{l, n, et). Now, we are in 
a position to explore the secular and periodic variations 
of Ca to the 3.5PN order, which will be done in the next 
two subsections. 



C. Secular variations 

First, let us extract the secular variations of Ca from 
Eqs. (30). This is achieved by using Eqs. (19a), which 
implies that the secular evolution of Ca can be obtained 
by orbital averaging the right-hand side of Eqs. (30), us- 
ing I. However, let us note the following points before 
we start. It is preferable to perform the orbital averag- 
ing in terms of u rather than /, as the right-hand sides 
of Eqs. (30) are explicit functions of u, and to this accu- 
racy we can use c?Z ~ (1 — et cosu)du, where we benefit 
from the fact that Eq. (27) for I does not have a IPN 



contribution. Moreover, the required integration over u 
can be easily facilitated by the following definite integral, 
available in Ref. [47], 

1 /■^'^ du 
2^ Jo (1 - gf cosu)^+i 

^ (l-e?)(^+i)/2^^(7T=^)' ^^^^ 

where P/v is the Legendre polynomial. 

These statements are mainly useful to obtain the differ- 
ential equations for n and et, namely, dn/dl and det/dl. 
Further, it is straightforward to express these differential 
equations in terms of the original time variable t rather 
than I, by using dl = ndt. In this way, we obtain dn/dt 
and det/dt, which explicitly read 



'dt ^ 5(1-6^)7/2 1^^ + 292e? 4- 37e^| + _^__^|20368 - 14784r; + (219880 - 15960077)6? 

+ (197022 - 141708r7)6t 4- (11717 - 828877)6?} , (32a) 

f - ^ Jil-'^flA '''-^'^'^] - ^^^^^^{340968- 22870477+ (880632 - 65125277)6? 

+ (125361 - 93184?7)6^| , (32b) 



where n and Ct, on the right-hand side of these equa- 
tions, stand for n and et, respectively. We have checked 
that our above results are in excellent agreement with 
equivalent expressions, available in Refs. [48, 49], com- 
puted using balance arguments involving local radiation 
damping and far-zone fluxes. As these secular evolutions 
of n and 6t, namely, n and et, are crucial for the phas- 
ing, we additionally obtain in Appendix B, for the first 
time, using balance arguments, 2PN accurate expressions 
for dn/dt and det/dt in harmonic coordinates, providing 
0{c~^) corrections. 

Now, let us turn our attention to the secular variation 
of Q and Ca, namely, q and cx. The arguments, employed 
in Ref. [25], to show that = = G\, where Gi = Fi/n 
and G\ = Fx/n, respectively, at the 2.5PN order are also 
extendable to the 3.5PN order. Further, we note that the 
right-hand sides of Eqs. (30c) and (30d) are functions of 



the form sin-u x /(cosw) and {v — u) x /(cosu), respec- 
tively, and hence they are odd under u —u. Therefore, 
their average over d/ ~ (1 — et cosu)du also exactly van- 
ishes, leading to G; = = G\ to the 3.5PN order. This is 
also consistent with another line of reasoning, presented 
in Ref. [25], that involves the time-odd character of the 
perturbing force yA' , dci/dv'', and 802/ dv^ , respectively, 
appearing in Eqs. (15), ending up with the conclusion 
that dci/dt and dc\/dt are time odd. Summarizing, we 
find that there are no secular evolutions for both ci and 
c\ to the IPN order of radiation reaction: 

^ = 0; ci{t) = ciito) , (33a) 
^ = 0; exit) = cxito) ■ (33b) 
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D. Periodic variations 

To complete this study, we focus now our attention on 
the differential equations for n, e^, c/, and c\, which give 



orbital period oscillations to our dynamical variables at 
0{c~^) and ©(c^''). First, let us consider the differential 
equations for n and it. They are 



dn _ %e'^nri f ^ _ 32 49 - 9e^ _ 35(1 - e^) \ _ e^^nrjr ^ 2 ^ 07 4\ _ -360 + 1176r? 

dl- 5 \x' X' J 5(l-e?)V2iy^ + ^y^'^*+^^^*/ 35 1 

+ ^680 - 1170477 ^ r _ ^^^2 34356^ + (36 - 756'n)e^] + [1470 - 47880?7 - (350 - 17080?7)e,2] 4r 

X X X 

+ [13510 + 3178077 - (24220 + 3052077)ef + (10710 - 126077)6^1 \ - (27594 + 588077)(1 ~ e?)^ 

X X 

11760(1 -e2)3] ^y^nri 



X^ / 280(1 e|)9/2 {^^^^^ ~ 1478477 + (219880 - 15960077)e4 + (197022 - 14170877)e( 

(11717- 828877)e?| , (34a) 



det 8^-^/3 1 _ g2 w 3 49_ge2 35 1 - 62] ^^/s £^^^^7 f r 

— = — - + ^ ^ , } + — ^ -^K^l 304+ 1216H + ^ - 4320 

dl 15et \x^ X^ X" X^ J 15(1-62)5/2 1 */ 3156* \ ^ 

+ 663677 + (5328 - 76Urj)e'f] ^ + [20430 - 3617677 - (25470 - 41216?7)6?] ^ + [ - 73650 + 106120?7 

+ (112692 - 116200?7)6? - (39042 - 10080?7)6^] ^ + [102312 - 15428077 - (201264 - 21616077)e? 

X^ 

+ (98952 - 618807?)6t] ^ + [2730 + 9534077 + (210 - 18690077)6^ - (8610 - 8778077)6^ + (5670 + 378077)e^] 

X X 

(82782 + 1764077)(1- 6,2)3 35280(1 - e^)^ ] ^ f/y, ^^^^^^^ 
7 + ? 1 + 2520(1 - 62)^2 (340968 - 228704,? 

+ (880632 - 651252?7)6? + (125361 - 93184?7)e^| , (34b) 

where n and e*, on the right-hand side of these equations, again stand for n and et, and we recall that the right-hand 
sides of Eqs. (34) are zero-average oscillatory functions of I. We have already argued that G; = = Gx to the IPN 
reactive order. This implies that the differential equations for e; and c\ are identical to those for q and c\, as given 
by Eqs. (30c) and (30d), but with n and et replaced by n and et, respectively. Symbolically, this reads 

dcj 

= RHS of Eq. (30c) [77 ^ n, 6t -> e*] , (35a) 
^ = RHS of Eq. (30d) [?7. ^ n, a Et] . (35b) 

One can analytically integrate Eqs. (34) and (35) to get n, et, cj, and eA as zero-average oscillatory functions of I. 
We find, when expressed in terms of it, 

^5/3717764 sin uf 602 6736? 314 - 2036? - llle^ 98- 12462 - 466^-^ 726^ 210(1-6^)3 
15(1-62)3 \ X ^ ? ^ ? ^ ? 

' ^'''"'^ r96 + 2926?H-376n(2tan-i('^^^#^)+6,sin. 



5(l-e2)V2l * n\ \1-I3tc0su, 

- ^ ' " ( [827796 - 6OI72O77 4- (4322828 - 313166077)6? -h (1584181 - 1132320?7)6fl - 
4200(1 — ei) I X 

- \pfniiQ - 37996077 + (1024628 - 73766077)6? - (1371149 - 99330077)e^ - (175755 - 12432077)e^] ^ 

X 

- [391116 - 33964077 - (296974 - 44072077)6? - (557800 - 66920?7)6^ + (442058 - 9744077)e? 

-(21600 - 705607?)6?] ^ -t- [196476 -f 155 4 007? - (401058 -f 89964077)6? 4- (24318 -f 1766520?7)6^ 
X 
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-(368634 - 145572077)e^ - (188370 - 433440?7)e?] — + [170856 - 504000ry - (668304 - 198576077)et 



-(964656 - 2903040?7)e^ - (592704 - 1834560r;)e^ + (110376 - 383040?7)e^ 4- (15 1 20 - 30 2 407/)et°] ^ 



_^ (115080 -I- 11760077)(l-e?)^ 201600(1 - e?)^ ] , C'^ 



nrj 



^6 ^7 J ' 280(1 - e?)9/2 {2°^^^ - ^^^^^'^ + (2^^^^° 159600^)e? 

(197022- 141708r;)e^-f (11717 - 8288r;)e^| (2 tan"i (tAt""") + ^tsmuX , (36a) 

J \ J- ~ Pi COS U J j 

^5/3,ysinu ri34+ 1069e?-h72e^ 134 -h 157e? - 291e^ 98 - 124e? - 46e^ -I- 72e? 210(1 - e?)^ 



45(l-e2)2\ X X* 

-^^^!^(304+ 121en(2tan-i ( ^f"" )+e,sinu] [78768 + 196 

15(l-e?)5/2l *i\ Vl-Acosuy J 37800(1 -e?)3\L 

(9997134 - 703346077)e? + (9942753 - 7434560?7)ef + (185760 - 13104077)e?] - 

X 

[78768 + 1960?7 + (4882614 - 360290077)et - (3266727 - 233422077)et - (1694655 - 1266720?7)e^] 
[337968 - 396200r/ + (1267458 - 233520r?)e? - (3136260 - 1131480r;)e^ + (1118274 + 22400?/) 



-f (412560 - 5241607?)e?] ^ + [-306432 + 7854007? + (2380266 - 3975720?7)e? - (5302206 - 7214760r7)e| 
-|-(4689342 - 5643960r/)ef - (1460970 - 16195207/)e^] ^ + [1419768 - 1512000?? - (6540912 - 5957280??)e^ 

-h(11965968 - 8709120??)e^ - (10850112 - 5503680?7)e^ + (4867128 - 114912077)e^ - (861840 + 9072077)6^^°] ^ 

(345240 -H352800?,)(l-e?)^ _ 604800(1 - e^)^ | ^je, r 

X^ X' } 2520(1 -e2)7/2 I ' 

-h (880632 - 651252?7)e2-f (125361 - 93184?;)e^||2tan-i (l^y^^) +etsinu|, (36b) 

2f/3?y/l44e^ , 18 - 258e? -56 + 92e^ - 36e^ 105(1 - e^)^ I . , 



456? I X X^ X^ 2(l-e?)3/2 

I (62640 - 141120?;)e^ ^ ^^^^ _ ^^^^^^^ _ ^^^^^^ _ gggg^Q^^^an J_ + T _ 80550 + 147700?] 
4725e|^ I X X 



(119310 - 6412007;)e2 - (21960 - 31500?7)e^] — + [215775 - 287175?] - (387450 - 599550?7)e? 



X 

-t-(17 1 675 - 3 1 23 757?)e^] \ + [-134316 + 233100?/ + (179802 - 454860?7)et2 + (43344 210420?7)e^ 
-(88830 - 11340,,,.?] - ("^7.5 . 44100„)(1 - e^P ^ T^^OOO^ ^ _^ j^^^^^ ^ 

-h (3486804 - 200018077)e? -I- (2246493 -I- 292180?;)e^ - (335790 - 968940?;)e^] I , (36c) 



45e 



2 



18 56- 36e? 105(1 - e?) 



^2 ^3 ^4 



134 + 103e? - 252e^) a/ 1 - - 134 - 295e? - 36e^ 



- 144e2 18 - 258et^ _^ 56 - 926^^ 36e4 105(1 - 



2\2 



^2 ^3 ^4 



4725e?(l - e?) 



21060 - 4977077 



2(l-e?)2 

(5400 + 4599077)e?] + \ - 45270 + 147700?? - (40290 + 1523207?)e? -I- (5760 -I- 8820?7)ef] 4t 
X X'^ 

[149625 - 287175?7 - (289800 - 438900?7)e2 + (140175 - 151725?7)e41 ^ + [ - 134316 4- 233100?? 

X 

(164682 - 454860?7)e? -t- (73584 + 210420??)e4 - (103950 - 11340?y)e6] - (118^55 + 44100??) (1 - eg)^ 



X X 
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75600(1 - ef )4 



[(73440 + 14112077)e? + (62640 - 14112077)et] - + [ - 21060 + 497707? 

X 



(64980 + 586740?7)e? - (178560 - 536970?7)e^] — + [45270 - 147700?? - (141900 - 78890077)6(2 
(118590 - 67270077)e^ - (21960 - 31500?7)e?] ^ + [ - 149625 + 287175?? + (470925 - 88672577)e 



X 



1 



(492975 - 911925?7)e^ + (171675 - 312375?7)e?] — r + [134316 - 233100?? - (314118 - 687960??)e? 

X 

(136458 - 665280?7)e^ + (132174 + 199080??)e? - (88830 - 11340?;)e?] ^ + (11^^55 + 44100?7)(1 - e^)^ 

X X 



2-\5 



1 



I [416448 + 1960?? + (3202044 - 200018077)e? + (1248573 + 292180?7)£ 



75600(1 -ef) 

8(l-e?)2 

(335790 - 968940?7)ef] Jl-ei - 416448 - I960?; - (4223964 - 1827140?7)e2 - (3117453 - 122402077)e 



(99810 -69300??)e^| 



(v — u) 



(36d) 



where (3t = (1 ^ ■\/l — ef)/et. The constant contribu- 
tions to the time evolution of c/ and c\, appearing in 
Eqs. (36c) and (36d), respectively, are required to guar- 
antee the zero-average behaviour. We note that the re- 
maining integral in Eq. (36d) can be numerically evalu- 
ated. 

The above results also modify the temporal evolution 
of the basic angles I and A, entering the reactive dynam- 
ics, Eqs. (5). According to Ref. [25], we see from the def- 
initions of l{t) and X{t), given by Eqs. (14), that we can 
also split these angles in secular and oscillatory pieces, 
denoted by /, A, and I, X, respectively, as 



where 



X{t) = X{t) + \[l;Ca{t)]. 



l{t) = I n{t)dt + ci{t) 



(37a) 
(37b) 



(38a) 



[1 + k{t)]n{t)dt + cx{t) . 



(38b) 



We note that ci(t) = c/(io) and c\{t) = c\{tQ) are con- 
stants [see Eqs. (33)]. The oscillatory contributions to I 
and A are given by 



' n 



n -n 
- + k- 
n n 



dl + cx{l), 



(39a) 
(39b) 



where k = (dk / dn)n + {dk/det)et denotes the oscillatory 
piece in k. 

Finally, to complete our study of the oscillatory contri- 
butions associated with the reactive dynamics, we com- 
pute the integrals of Eqs. (39) and add them to the pre- 
vious results for Ciil) and C\{1), respectively. We find 



IQ'] Co) 



(602 + 673e^)x + (314 - 2036^ - lllef ) In x - (602 + 673e^) 



-98 + 124e2 + 46e4 - 72e 



5(1 - eiyn 



96 + 292ef + 37e; 



2 tan 



Pt sin u 
1 — j3t cos u 



et sm u 



Xdu 



15(1 -e2)3 
105(1 - e2)3 
? 

^ 4200(1 - e2)4 

+ [522276- 379960?y+ (1024628- 7376607?)e^ - (1371149 - 99 3 3 007?)e^ - (175755 - 124320?7)e^] Inx 
- 827796 + 601720?; - (4322828 - 313166077)e? - (1584181 - 1132320??)e* + [ - 391116 + 339640?? 

+ (296974 - 440720?7)e? + (557800 - 66920?7)ef - (442058 - 97440?7)e^ - (21600 - 7056077)e?] - 



[827796 - 601720?? + (4322828 - 3131660?7)e^ + (1584181 - 11323207?)e^] x 



98238 - 77700?; + (200529 + 449820?;)ej2 - (12159 + 883260?7)e^ - (184317 - 727860?;)e^ 
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(94185 - 21672077)6^1 + \ - 56952 + 168000ry + (222768 - 661920r;)e,2 - (321552 - 967680r7)e 

X 



(197568 - 611520r7)e^ - (36792 - 127680?7)e^ - (5040 - 1008077)e 



101 



(28770 + 2940077)(1 - 6(2) 



2\5 



40320(1 - e2 



A(Z;ca) 



+ (11717- 828877)e^|| J 



280(1 - e2)9/2 

1 f Ptsmu 



1 20368 - 14784r; + (219880 - 159600r7)e? + (197022 - 141708?7)e: 



2 tan 



1 — (it cos u 



et sill u 



Xdu} + ci{l) , 



(40a) 



15(1 



(602 + 673ei)x + (314 - 2036^^ - llle^) In x - (602 + 673e^) 



-98 + 124e? + 46ef - 72e? 



105(1-6^)3 



5(l-e?)V2 



96 + 292e? + 37ef 



2 tan 



Pt sin u 
1 — /3t cos u 



X 



et sm u 



Xdu 



[1595556 - 601720?7 + (4666388 - 313166077)e? + (1543861 - 1132320?7)e^] x 



4200(1 - el, 

+ [886836 - 37996077 + (652508 - 737660?7)e2 - (1363589 - 99330077)e^ - (175755 - 12432077)e?] Inx 
- 1595556 + 60172077 - (4666388 - 3131660?7)e2 - (1543861 - 113232077)e^ + [ - 473436 + 339640?7 

+ (401134 - 440720?7)e? + (596440 - 66920?7)ef - (502538 - 97440?7)e? - (21600 - 7056077)e?] - 

X 

+ [ - 186438 - 7770077 + (465129 + 44982077)6? - (276759 + 88326077)e^ - (96117 - 7278607;)e^ 

+ (94185 - 21672077)6^1 \+\- 56952 + I68OOO77 + (222768 - 66192077)6^ - (321552 - 9676807;)ef 

X 

+ (197568 - 6115207;)e^ - (36792 - 1276807;)e^ - (5040 - 1008077)64^°] ^ - (28770 + 2940077)(1 - ) 

X X 



40320(1 



' ' tttttttt-^-^Wt^ (47248 - 1478477 + (267592 - 1596007;)et + (193830 - 14170877)6^ 
280(1 — 6( )9/2 L 



+ 



(11717- 828877)6^1 



2 tan 



I3t sin u 
1 — /3t cos u 



+ 6t sm u 



Xdu\ + cx{l) 



(40b) 



r 



where ci{l) and c\{l) are given by Eqs. (36c) and (36d), 
respectively, and Pt = (1 — ■\/l — ef)/et. Note that the 
contributions to A(/) arising from the periastron advance 
constant k only appear at 0{c~''). 

In the next section, we plot the current results obtained 
in these subsections and their influences on /i+ and /ix- 



IV. VISUALIZATION OF SOME EXEMPLARY 
3.5PN ACCURATE RESULTS 

In this section, we present a few samples of the tem- 
poral evolution of Cq, and Cq, . These "varying constants" 
cause secular and periodic variations in the dynamical 
variables that appear in the expressions for /7+(t) and 
hx{t), as given by Eqs. (3). 

When we evolve orbital elements and gravitational 
waveforms, we have to make sure that the eccentric or- 
bits we study lie inside the area of validity of our ap- 
proach. This is guaranteed by terminating the orbital 
evolution when the 3PN accurate version of the following 
constraint, discussed in Ref. [25], connecting the param- 



eters 77 (via ^ = GMn/c') and 64, is violated: 



1 



(1 - e2)3/2 



< 



1 



483/2 



3.007 X 10- 



(41) 



where the dimenssionless angular-momentum variable is 
given by j = cL/{^GM). The above inequality, obtain- 
able from Eqs. (20) and (21) in Ref. [25], ensures that 
we are staying sufliciently far away from the last stable 
orbit (LSO). In this way, the usage of Eqs. (7) and (8) for 
the conservative orbital dynamics and their counterparts 
when the dynamics is reactive is fully justifled. Accord- 
ing to Eq. (41), we will terminate the orbital evolution 
when j = V48. [Recall that for a test particle around a 
Schwarzschild black hole, the LSO is at j = vT2.] Go- 
ing beyond the above restriction will require the effective 
one body approach, being developed by Damour and his 
collaborators [see related discussions in Ref. [25]]. 

We plot in this section only dimensionless quantities 
like ^, in terms of dimensionless variables by making use 
of the following inherent scaling. The conversion to famil- 
iar quantities like orbital frequency / (in hertz) is given 
by / = 77/(27r) = c^^/{2ttGM) = 3.2312 x 10^^(Mq/M). 
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Figure 1: The plots for n/rii and h/n versus 1/(2-k), which 
gives the number of orbital revolutions. The adiabatic in- 
crease of n is clearly visible in panel 1, and the quasi-periodic 
nature of the variations in n is portrayed in panels 2-6. These 
variations are governed by the reactive 2.5PN and 3.5PN 
equations of motion. In the second and third ro"w, these con- 
tributions to h are plotted individually and separated for the 
initial and final stages. The parameters ej and e{ denote ini- 
tial and final values of the time eccentricity Ct, "while and 
stand for similar values of the adimensional mean motion 
^ = GMn/(? . The panels are plotted for rj = 0.25 and the 
orbital evolution is terminated when j = 






10 



Figure 2: The plots for the scaled h+(t) and hx (t) (Newtonian 
in amplitude and 3.5PN in orbital motion) as functions of 
1/{2-k). The slow chirping and the amplitude modulation due 
to the periastron precession are clearly visible in the two upper 
panels. In the two bottom panels, we zoom into the initial 
stages of the orbital evolution in order to show the effect of 
the periodic orbital motion and the periastron advance on 
the scaled h+(t) and hx{t). The initial and final values of the 
relevant orbital elements are marked on top of the plots. The 
panels are plotted for a binary consisting of equal masses, so 
that 7] = 0.25, and the orbital inclination angle is given by 
i = tt/S. The orbital evolution is terminated when j -- 



This implies that for a compact binary with the total 
mass M = Mq and ^ = 10^"^, the orbital frequency "will 
be - 30 Hz. 

In Fig. 1 we plot n/rii, -where is the initial value 
of n, and h/n, as functions of 1/(2^), -which gives the 
evolution in terms of elapsed orbital cycles. We clearly 
see an adiabatic increase of h as -well as the quasi-periodic 
variations of h. These variations are governed by the 
reactive 2.5PN and 3.5PN equations of motion. In the 
second and third ro-w of Fig. 1, these contributions to n 
are plotted individually and separated for the initial and 
final stages. These 2.5PN and 3.5PN contributions are 
obviously in-phase, and -we observe that the scaled 3.5PN 
contributions are only by a factor of ^ 10 smaller than 
their scaled 2.5PN counterparts. 

Though, -we have all the required computations to plot 
the secular and quasi-periodic variations in ej, q, and cx 
to the IPN reactive order, we do not attempt it here. As 
expected, we observed the same features at the 3.5PN 
order, as detailed in Figs. 2 and 3 in Ref. [25] at the 
2.5PN order — the adiabatic decrease of et, the periodic 
variations of it, no secular evolution of c/ and c\, but 
periodic variations in q and cx — and that is the main 
reason for not duplicating these figures. 

Finally, we plot in Fig. 2 the scaled h^{t) and hx{t), 
evolving under gravitational radiation reaction, as func- 



tions of V(27r). [We factored out Gn/{c^R') appearmg 
in H+Iq and hxlq to get the scaled waveforms.] We em- 
ploy for these figures polarization amplitudes, which are 
Newtonian accurate, as given by Eqs. (3), while the or- 
bital motion is 3.5PN accurate. We clearly see "chirping" 
due to radiation damping, amplitude modulation due to 
periastron precession, and also orbital period variations. 

We note that Figs. 1 and 2 can be used to illustrate the 
various aspects of a compact binary inspiral from sources 
relevant for both LIGO and LISA. This is based on the 
above mentioned scaling argument. Let us detail this in 
case of the following two scenarios. For instance, if we 
choose for the total mass M ~ 2.8Mq — a binary inspi- 
ral involving two IAMq compact objects — the variation 
of ^ from 2.069 x lO'^ to 3.0209 x lO^^ in - 264 orbital 
cycles corresponds to a increasing orbital frequency from 
150 Hz to 219 Hz in ~ 7.58 s. Similarly, for the choice of 
M ~ IO^Mq — a binary inspiral involving two supermas- 
sive black holes — the variation of ^ from 2.069 x 10""^ 
to 3.0209 X 10"'^ in ^ 264 orbital cycles corresponds to 
a increasing orbital frequency from ~ 4.2 x 10"'^ Hz to 
- 6.1 X 10-3 Hz in - 3.1 days. 

We conclude by noting that Fig. 1 shows clearly the 
existence of periodic variations in the orbital elements — 
analytically investigated for the first time in Ref. [25] and 
improved to the 3.5PN order in this paper. 
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V. CONCLUSIONS 

Let us recapitulate. In this paper, we have incorpo- 
rated the 3PN accurate conservative and the IPN accu- 
rate reactive dynamics in harmonic coordinates into the 
phasing formaHsm to the 3.5PN order as a natural exten- 
sion of the work presented in Ref. [25]. This extension 
was possible due to the very recent determination of the 
3PN accurate generalized quasi-Keplerian parametriza- 
tion for the conservative orbital motion of nonspinning 
compact binaries in eccentric orbits [31]. We applied 
the method of Ref. [25] to construct, almost analytically, 
templates for GW signals emitted by compact binaries 
moving in inspiralling and slowly precessing eccentric or- 
bits. An improved method of variation of arbitrary con- 
stants, explained in great detail in Ref. [25], allowed us 
to combine the three different but relevant time scales, 
namely, those associated with the radial motion (orbital 
period), advance of periastron, and radiation reaction, 
to the 3.5PN order without making the usual approxi- 
mation of treating adiabatically the radiative time scale. 
In this context, we recall that the two-scale decompo- 
sition helped us to model accurately and efficiently the 
time evolution of the associated dynamical variables. We 
employed harmonic coordinates in this paper as calcula- 
tions that provided search templates for compact bina- 
ries in quasi-circular orbits usually employ the harmonic 
gauge. 

The expHcit computations provided in this paper will 
be required to construct accurate and efficient search 
templates for gravitational waves from compact binaries 
of arbitrary mass ratio moving in inspiralling eccentric or- 
bits. Our detailed calculations will have to be employed, 
if the earth-based GW interferometers plan to search for 
gravitational waves from compact binaries with residual 
eccentricities, motivated by a plethora of recent astro- 
physical investigations [4-13]. The proposed space-based 
GW interferometers Hke LISA, BBO, and DECIGO will 
have to depend on our results to do astrophysics. It is 
interesting to note that our results will be required by 
LISA to search for gravitational waves from stellar-mass, 
intermediate-mass, and supermassive black-hole binaries 
as these binaries will Hkely to be in inspiralling eccen- 
tric orbits. Another area where our computations can be 
quite effective will be the early stages of extreme mass 
ratio inspiral (EMRI) as relevant for LISA. Our current 
results should be also useful to benchmark efforts that 
are required to obtain reliable EMRI templates [50]. 

There are many avenues that will require detailed in- 
vestigations in the near future and we list only a few of 
them below. In this paper, the conservative dynamics 
was restricted to compact binaries consisting of nonspin- 
ning point masses. Naturally, it is desirable to include 
spin effects into our computations. A first step in this 
direction was taken in Ref. [43], where the 3PN accurate 
generalized quasi-Keplerian parametrization for the con- 
servative dynamics of spinning compact binaries, moving 
in eccentric orbits, when the spin effects are restricted 



to the leading-order spin-orbit interaction, is presented. 
We also neglected, for simplicity, explicit PN corrections 
to the GW polarization amplitudes of /i+ and h X , and 
restricted them to their leading quadrupolar order. How- 
ever, it is possible to obtain 2PN accurate corrections to 
these amplitudes, using Refs. [25, 41]. In order to make 
the numerical implementation of our computations more 
efficient and accurate, it is also desirable to provide better 
ways of solving the 3PN accurate Kepler equation. An- 
other line of investigation should deal with a extension of 
these computations so that we have a dependable descrip- 
tion for the orbital evolution near the LSO. Finally, these 
templates naturally trigger lots of data analysis investi- 
gations relevant for both ground-based and space-based 
GW interferometers. Many of the above mentioned issues 
are currently under investigation. 
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Appendix A: CONSTRUCTION OF AN EXACT 
RELATION FOR v-u 

In this appendix, we provide the details involved in the 
derivation of the exact relation for v — u, which is also 
periodic in u: 



V — u = 2 tan 



, smu 



(Al) 



where f3^ = {1 — y 1 — e^)/e0. This relation allows us to 

avoid the usage of the commonly employed infinite series 
expression for v — u, namely. 



V — u ^ 2 sm lu . 



(A2) 



In order to deduce Eq. (Al), we start from the follow- 
ing identity 



w — u 2 tan ^ 
Using therein 

tan(a — /3) 



tan 



tan a — tan (3 
1 -|- tan a tan (3 ' 



(A3) 



(A4) 
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leads to 



w — u = 2 tan ^ 



tan I — tan ^ 



1 + tan I tan | 



(A5) 



The relation connecting the true anomaly v to the eccen- 
tric anomaly u, as given by Eq. (9c), is used to replace 
tan(w/2) in the above equation. In this way, we obtain 



V — u = 2 tan 



(g0-l)tanf 



I + Q4, tan- 



2 u 
2 



where 



1 + 

1- Cd 



1/2 



With the help of 



u smu 

tan — = 

2 1 + cos u 



we rewrite Eq. (A6) as 

V — u ~ 2 tan^"'^ 
Now, let us call 



smu 



6-1 



cosu 



(A6) 



(A7) 



(A8) 



(A9) 



heuristic arguments, detailed in Ref. [48] and in Sec. VI 
in Ref. [25]. In the heuristic determination of the evolu- 
tion equations for n and et , one employs PN accurate ex- 
pressions for n and et, and the far-zone (FZ) energy and 
angular-momentum fluxes. The PN accurate expressions 
for dn/dt and det/dt are then obtained by differentiating 
the PN accurate expressions for n and et, expressed in 
terms of E and L, with respect to time and then heuris- 
tically equating the resulting time derivatives of E and 
L to the orbital averaged expressions for the FZ energy 
and angular-momentum fluxes. For the ease of imple- 
mentation, we split the 2PN accurate computations of 
dn/dt and det/dt into two parts. The flrst part contains 
the purely "instantaneous" 2PN corrections and the sec- 
ond part considers the so-called "tail" contributions [51], 
appearing at the 1.5PN (reactive) order and derived for 
the flrst time in Refs. [52, 53]. The computations to get 
the instantaneous contributions begin with the 2PN cor- 
rections to the FZ fluxes, in harmonic gauge, in terms 
of r, f, and available in Ref. [41]. These FZ fluxes 
are orbital averaged, using the 2PN accurate general- 
ized quasi-Keplerian parametrization for elliptical orbits 
in harmonic gauge, following the prescripton detailed in 
Ref. [48]. We perform the orbital average by using an 
additional ingredient, namely, the relation connecting dl 
and du to 2PN order in harmonic coordinates 



1 ' 



which can be simplifled to 



(AlO) 



(All) 



Finally, the combination of Eqs. (A9)-(A11) directly 
leads to Eq. (Al). 

In addition, analog to the above derivation, we con- 
structed the following exact relation, involving the time 
eccentricity et instead of e^. 



2 tan 



tan — 

1-et/ 2 



2 tan" 



Pt sin u 
1 — /3t cos u 



(A12) 



where /3t = (1 — ^1 — )/et, which was required in case 
of Eqs. (36a) and (36b). 



dl_ 
du 



^4/3 



p2 



{15v-v'){l-e'tf^' 



+ (60 - 39?7 + ri^)x\ll-el - (60 - 2477)x^ 



(Bl) 



where x = 1 ^ cos u. The resulting deflnite integrals 
are easily computed, using Eq. (31). Now, we compute 
the time derivatives of the PN accurate expressions for n 
and et and equate the resulting time derivatives of E and 
L to the orbital averaged expressions for the FZ energy 
and angular-momentum fluxes, respectively, to get the 
PN accurate expressions for dn/dt and det/dt in terms of 
E, L, M, and r]. Finally, we use Eqs. (21) to obtain the 
differential equations for fi and et in terms of fi, et, M, 
and 77. The resulting 2PN accurate instantaneous con- 
tributions to dn/dt and det/dt, in harmonic coordinates, 
are given by 



Appendix B: 2PN ACCURATE ADIABATIC 
EVOLUTION OF n AND et IN HARMONIC 
COORDINATES 

Following Ref. [25] , let us briefly show in this appendix 
how to obtain, in harmonic coordinates, the 2PN accu- 
rate secular changes in n and et. The PN accurate dif- 
ferential equations for n and et are computed using the 



^ = f {^N ^ ^IPN _^ ^2PN| ^ (g2a) 

^ = --e/^n^-et {et + e^™ + e1™} , (B2b) 

where the various instantaneous PN accurate corrections, 
namely, n^, n^^^, n?^^ , , e\^^ , and ef^^ , read 
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^IPN 



ilPN 



5(l_e-,2)7/2 {96 + 292e1 + 37el} , (B3a) 



t , 

"2/3 



280(1 e?)9/2 {^^^^^ ~ 14784r; + (219880 - \hm^Q^)e\ + (197022 - 141708?7)ef + (11717 - 8288r;)e^| , 

(B3b) 

2PN = / ^ 1 12592864 - 13677408r; + 19031047;^ + (133049696 - 1855385287? + 612820327;2)e2 

30240(1 — e^Y^I ^ y 

+ (284496744- 411892776/;+ 16650606077^)6^ + (112598442 - 14208906677+ 6482884877^)6^ 



f (3523113 - 325998077 + 196425677^)6? + 3024(96 + 42686^ + 4386e^ + 175e?)(5 - e?| , (B3c) 

^^^^^{304+ 1216-?}, (B3d) 

^^^^^^ — rjy^ |340968 - 22870477 + (880632 - 651252?7)e? + (125361 - 93184?7)e^| , (B3e) 

*™ = 30240(^ ' -2)9/2 (20815216 - 25375248?? + 454809677^ + (87568332 - 128909916?? + 48711348?7^)e^ 
+ (69916862 - 9352257077 + 428100967?2)et + (3786543 - 4344852?? + 2758560??^)e^ 



+ 1008(2672 + 6963e? + 565e^)(5 - 2??)y^l -e?| , (B3f) 



I 

where ^ = GMfi/c^ . The tail contributions to dfi/dt and curate contributions to dn/dt and det/dt in ADM gauge. 

det/dt, which appear at the 1.5PN order, are already Note also the expected differences at higher-PN orders 

presented in Sec. VI in Ref. [25], given by Eqs. (70) and between the harmonic and the ADM gauge. We con- 

(71) therein. elude by noting that the article providing the 3PN accu- 

We have checked that to the IPN order the above con- rate contributions to dn/dt and det/dt is currently under 

tributions are in excellent agreement with Eqs. (68) and preparation [54]. 
(69) in Ref. [25], which give the instantaneous 2PN ac- 
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